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We discuss the gamma-ray signal from dark matter annihilation in our Galaxy and in external 
objects, namely the Large Magellanic Cloud, the Andromeda Galaxy (M31) and M87. We derive 
predictions for the fluxes in a low energy realization of the Minimal Supersymmetric Standard Model 
and compare them with current data from EGRET, CANGAROO-II and HEGRA and with the 
capabilities of new-generation satellite-borne experiments, like GLAST, and ground-based Cerenkov 
telescopes, like VERITAS. We find fluxes below the level required to explain the possible indications 
of a 7-ray excess shown by CANGAROO-II (toward the Galactic Center) and HEGRA (from M87). 
As far as future experiments are concerned, we show that only the signal from the galactic center 
could be accessible to both satellite-borne experiments and to ACTs, even though this requires very 
steep dark matter density profiles. 
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I. INTRODUCTION 

The nature of the Cold Dark Matter (CDM) which is 
believed to compose galactic halos is probably the most 
important open issue in present Cosmology. A popu- 
lar solution to this puzzle is given by the lightest super- 
symmetric particle (LSP) which, in most supcrsymmetry 
breaking scenarios, is the neutralino x- I n this case Dark 
Matter (DM) would be not so dark after all, since x~X 
annihilation is expected to lead, among other final states, 
to a 7 signal which could in principle be detected above 
known backgrounds. In particular, since the neutralino 
annihilation rate is proportional to the square of its den- 
sity, a signal enhancement is expected in high density 
regions like the center of our Galaxy or that of exter- 
nal ones, with the exciting possibility that such 7-rays 
might be identified by forthcoming or just operating at- 
mospheric Cerenkov telescopes (ACT) such as VERITAS 
Q HESS and MAGIC Q or by satellite-borne detec- 
tors like GLAST let alone the even more intriguing 
chance that a hint of an exotic source of 7-rays could 
actually be already present in the data of existing exper- 
iments, like EGRET Q or CANGAROO-II 0. However, 
assessing the size of such signals depends on many un- 
certain aspects of both astrophysics and particle physics. 
For instance, the central structure of the DM halos is far 
from being well determined, and this can lead to uncer- 



tainties in the calculation of expected 7 rates spanning 
several orders of magnitude. Another sensitive issue is 
the presence of substructures in galactic halos, which can 
change predictions as compared to a smooth mass distri- 
bution. 

The aim of the present paper is to investigate the pos- 
sibility that neutralino annihilations in the halo of our 
galaxy @, H S E3, or that of external ones [U El 
(namely the Large Magellanic Cloud, the Andromeda 
Galaxy and M87) could produce detectable fluxes of 7- 
rays. To this purpose we will discuss present astrophysi- 
cal uncertainties and focus on deriving consistent predic- 
tions for these fluxes in a specific realization of supcrsym- 
metry, the effective Minimal Supersymmetric Standard 
Model (MSSM). 

The plan of the paper is as follows: in Section II the 
main ingredients for the calculation of the 7-ray flux 
from neutralino annihilation are introduced; in Section 
III we discuss the contribution to the flux calculation 
coming from astrophysics, while in Section IV the contri- 
bution from particle physics is discussed, and the effective 
MSSM Supersymmetric model is outlined. In Section V 
we show our results and compare them to present data 
and the prospects of future experiments; finally, Section 
VI is devoted to our conclusions. 

II. THE 7-RAY FLUX 



The diffuse photon flux from neutralino annihilation in 
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and observed by a detector with angular resolution 6 can 
be written as: 
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The energy dependence in Eq. is given by the anni- 
hilation spectrum: 
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where (<T am if) is the neutralino self annihilation cross- 
section times the relative velocity of the two annihilating 
particles, dN? /dE 7 is the differential photon spectrum 
for a given /-labeled annihilation final state with branch- 
ing ratio Bf and m x denotes the neutralino mass. The 
geometry-dependence is given by the line-of-sight inte- 
gral, defined as: 
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for the diffuse emission of our Galaxy, and 
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for the emission from an extragalactic object located at 
the direction i/j. In Eq. @, p x (r) is the dark matter den- 
sity profile, r is the galactocentric distance, related to the 

distance A from us by r 



1 X 2 + Rq - 2Ai?Q cosV (Rq 
is the distance of the Sun from the galactic center) and 
AQ(ip, 9) is the solid angle of observation pointing in the 
direction of observation ip and for an angular resolution 
of the detector 8. Moreover, in Eq. |J3J d is the distance 
of the external object from us, Rg is the radius of the ex- 
ternal galaxy and r max (Af2) is the maximal distance from 
the center of the external galaxy which is seen within the 
solid angle Afl(ip,6). 

We focus our attention on the fact that Eq. Q is fac- 
torized into two distinct terms: a "cosmological factor" 
<j)Cosmo w j 1 j c j 1 takes into account the geometrical distri- 
bution of DM in the Universe, and a "supersymmetric 
factor" <J> SUSY which contains the information about the 
nature of dark matter. In Sections IIIII and IIVI we will 
present results on the two factors separately. 



III. THE "COSMOLOGICAL FACTOR" 

In the following we present the determination of the 
"cosmological factor" $ cosm ° ; a s defined in Eqs. © and 
(0J. The dependence of <j> cosmo on the astrophysical and 



cosmological details that we explore here is based on the 
determination of the shape of the dark matter halo. This 
takes into account the possible existence and prominence 
of central cusps, the study of the physical extent of the 
constant-density inner core, and the possible presence of 
a population of sub-halos. We remind the reader that, 
for the moment, no definitive answer can be given to 
these questions by experimental constraints. In particu- 
lar, the discussion about the possible existence of a halo 
with a cuspy behavior in its inner regions is still quite 
open. Moreover, theoretical predictions differ substan- 
tially among themselves, or take into account different 
input parameters. 

These facts reflect themselves in a large uncertainty 
in the predictions of the gamma-ray fluxes arising from 
$ cosmo , as it is discussed and quantified in the following. 

A. Modeling the Dark Matter Halo 

The modeling of the DM density profile is an open 
question. It can be addressed through numerical N- 
body simulations whose scale resolution is about few 
xl0 _3 rioo, where rioo is defined as the radius within 
which the halo average density is about 100p c (p c is the 
critical density). The very inner slope of the profile is 
then usually just extrapolated and does not take into ac- 
count interactions with the baryons which fall in the DM 
potential well. A number of profiles have been proposed. 
Here we discuss some of the profiles which arc compatible 
with observations and which we will use in our analysis. 

In our calculation we mainly focus on the NFW profile 
(hereafter NFW97) Q 
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and the Moore et al. profile (M99) Q : 
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The scale radii r\ and the scale densities p\ (i = 
NFW97, M99) can be deduced by observations (the virial 
mass of the halo or the rotation curves) and by theoret- 
ical considerations that allow to determine the concen- 
tration parameter c = r V i r /r s (the virial radius r v j r is 
defined as the radius within which the halo average den- 
sity is 200p c ). The concentration parameters, cnfwq7 
and Cmoo = 0-64 cnfwq7> have been computed accord- 
ing to Ref. 0| with the assumption of a CDM power 
spectrum with a shape parameter T = 0.2 normalized to 
cr 8 = 0.9. 
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FIG. 1: Comparison between cuspy and cored dark matter 
density profiles for the Milky Way, as a function of the dis- 
tance from the center of the Galaxy. All the curves are nor- 
malized to po = p(Rq) — 0.3 GeV cm -3 . 



In addition to the two profiles mentioned before, we 
include in our predictions the conservative modified 
isothermal profile with a constant density core (iso-core): 
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and a profile which has been recently proposed by Moore 
and collaborators (M04) Q: 
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Fig. ^ shows the comparison among the above men- 
tioned profiles for the Milky Way, normalized to a local 
density of 0.3 GeV cm -3 and to i? Q =8.5 kpc. Two more 
profiles are shown for comparison on Fig. ^ One is the 
numerical profile obtained in Ref. [T^ | when the adiabatic 
growth of a central black hole is taken into consideration 
(adiab-NFW). This hypothesis of black- hole formation 
has been applied here to the NFW97 profile, and the den- 
sity profile has been normalized as previously mentioned. 
The resulting profile has a behavior at the galactic center 
which is similar to the one of the M99 profile, therefore 
we won't discuss it in more details. The last profile which 
is shown in the figure is a cored one recently obtained in 
Rcf. Q (N03): 
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Galaxy 


mass {Mq) 


distance (Kpc) 


TVir (Kpc) 


MW 


1.0 ■ 10 12 


8.5 


205 


LMC 


1.4 ■ 10 10 


49 


49 


M31 


2.0 ■ 10 12 


770 


258 



TABLE I: Masses, distances and virial radii for the Milky 
Way, the LMC and M31. 



Profile 


scale radius r s 


(Kpc) 


scale density p s 


(M B kpc- 3 ) 


NFW97 


21.746 




5.376 ■ 


10 6 


M99 


34.52 




1.060 ■ 


10 6 


M04 


32.625 




2.541 ■ 


10 6 


iso-core 


4 




7.898 ■ 


10 6 



TABLE II: Scale radii and scale densities for the NFW97, 
M99, M04 and iso-core density profiles calculated for the 
Milky Way. 



where a = 0.17, rf ° 3 = rf FW97 and p™3 = / ,nfw97/ 4 _ 
As noticed in Ref. [lfij . this profile is compatible with 
the M04 as far as the resolution of the N-body simulation 
holds. In the inner part of the Galaxy, it is an extrapola- 
tion which postulates the existence of a constant density 
core. Another recently proposed profile which docs not 
exhibit singular behaviour, and which has been shown to 
be able to reproduce to a good precision the rotational 
velocities of low surface brightness galaxies ! is given 
in Ref. j^. Predictions of gamma-ray fluxes for this 
profile are given in Rcf. 

Since profiles shallower than the NFW97 hardly give 
observable fluxes of photons, we will not discuss it in 
detail. Studying the cored halos, we will limit ourselves 
to the iso-corc profile, which is pretty conservative. 

Integrating the squared density along the line of sight 
introduces divergences when cuspy profiles are consid- 
ered. Therefore we enforce a cut-off radius r cut to the 
density profile, with a constant density core therein. The 
smallest value for the cut-off radius which we will use is 
^cut = 10~ 8 kpc, a value we will discuss in the next Sec- 
tion, where the effect of varying r cut , both for our Galaxy 
and for the external ones, will be discussed. 

The analysis of Ref. [Llj shows that a number of ex- 
ternal galaxies shine above the Galactic foreground. In 
the following we will focus on the two most prominent 
galaxies at large angles with respect to the Galactic Cen- 
ter, namely the Large Magellanic Cloud (LMC) and the 
Andromeda Galaxy (M31) |n| Table d shows the as- 
trophysical parameters for the Milky Way, the LMC and 
M31, while Tables ITT1 ITTT1 and ITVl show the scale radius 
ans the scale density parameters used in our calculations. 
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Profile 


scale radius r s 


(Kpc) 


scale density p s 


(M kpc- 3 ) 


NFW97 


4.353 




8.50 • 


10 6 


M99 


6.8 




1.80 ■ 


10 6 


M04 


6.426 




3.22 • 


10 6 


iso-core 


1.5 




2.17 ■ 


10 7 



TABLE III: Scale radii and scale densities for the NFW97, 
M99, M04 and iso-core density profiles calculated for the 
LMC. 



Profile 


scale radius r 3 (Kpc) 


scale density p a (Mq kpc 3 ) 


NFW97 


30.271 


4.20 ■ 10 6 


M99 


47.298 


0.86 • 10 6 


M04 


44.697 


1.55 • 10 6 


iso-core 


4 


7.898 ■ 10 6 



TABLE IV: Scale radii and scale densities for the NFW97, 
M99, M04 and iso-core density profiles calculated for M31. 



1. Comment on the experimental constraints on the inner 
part of galaxies 

As we have seen, theoretical estimates of the inner 
slope a of the DM density profile p{r) oc r~ a are still un- 
certain. Moreover, observations which should constrain 
the a parameter do not give clear and definitive answers 
on its value. A number of works give in fact non-unique 
values for the slope. 

In Ref. jijj spatially resolved spectra of the diffuse 
hot (X-rays) gas of galaxies and clusters measured with 
the Chandra satellite were used to infer the radial mass 
distribution of the considered systems. An analysis was 
done on 2 clusters which are relaxed in their cores on 
O(10 2 kpc) to O(Mpc) scales and do not have strong 
radio sources in their center. Resulting values for a are 
1.25 and 1.35. A value of a less than 1 is found when 
disturbed X-ray surface brightness clusters are used. Yet 
the X-ray method uses the double assumption of a single 
phase gas in hydrostatic equilibrium, which for instance 
is questionable in the central regions where rapid cooling 
occurs. 

Other studies of radial mass profiles inferred by the ra- 
dial profile of the intracluster medium density and tem- 
perature measured with Chandra can be found in Ref. 
[22j where the analysis of 5 clusters give s 1 < a < 2. 

Different results are found by Ref. |23j using Low Sur- 
face Brightness (LSB) Galaxies rotation curves. Fits to 
their measured curves give a mean value (a) = 0.2, al- 
though tails in the distribution extend further, up to 
a — 2. In Ref. [3 a combination of strong-lensing data 



and spectroscopic measurements of stellar dynamics of 
the brightest cluster galaxies was used to derive values of 
a. Three clusters, containing both radial and tangential 
arcs, have been found. The obtained distribution gives 
(a) = 0.52 with Aa = 0.3. 

In Ref. ji^ the full radial extent of LSB galaxies ro- 
tation curves, instead of its inner portion, was used to 
determine the inner slope of the DM density profile. Con- 
vergence criteria for the N-body simulations taken from 
Ref. [2rj | give a minimum radius for which simulations 
are reliable r conv = l/i -1 kpc. It is shown that, at that 
radius, 2/3 of the sample in Ref. ji^l is consistent with a 
profile which lies between the simulated NFW97 and M99 
ones. There are inconsistencies with CDM predictions in 
those galaxies which show a sharp transition between the 
rising and flat part of the rotation curve. This is due to 
the fact that rotation curves of gas disks are compared 
with the spherically-averaged circular velocity profiles of 
DM halos. This assumption may not be correct in non 
regular galaxies. 

Another study of high resolution H a rotation curves 
for dwarf and LSB galaxies has been recently carried 
out in Ref. j^- In that work it is shown that rota- 
tion curves data are insufficient to rule out halos with 
a = 1, although none of the galaxies require an inner 
cuspy profile instead of a core density feature. Results 
on a range from to 1.2, although the quality of the 
fit is good only up to a = 1. Other analysis on large 
sets of data of high-resolution rotation curves also show 
consistency with cored mass distributions [3] . 

An indirect estimate of a can be inferred through the 
weak gravitational lensing measurements of X-ray lumi- 
nous clusters |29$ : one finds 0.9 < a < 1.6. 

The analysis of the microlensing optical depth toward 
the Galactic Center was performed in Ref. |3fJ . Assum- 
ing a naive spherically symmetric profile normalized to 
our position in the Milky Way, the authors find a = 0.4. 
They argue that the value a = 1 can be reached by con- 
sidering a flattened halo with a ratio of polar to equato- 
rial axis of 0.7. 



B. Including the Effect of the Inner Core 

There exists a physical minimal radius, r cut , within 
which the self-annihilation rate £; ~ {(<J aan v) nyf rmt))" 1 
equals the dynamical time td yn ~ [Gp)~^ [3jJ, where 
p is the mean halo density and n x is the ncutralino 
number density. When this procedure is applied to the 
density profiles we are using, the evaluated r cut are of 
the order of 10 -8 — 10~ 9 kpc for the M99 profile and of 
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1(T 13 - 1CT 14 kpc for a NFW97. Evaluating the constant 
core is indeed a much more complicate issue. Taking 
into account additional effects, especially tidal interac- 
tions, the central core of galaxies can significantly ex- 
ceed the values quoted above, reaching values as large as 
0(0.1 — l)kpc |3J|. We want to remind that also nu- 
merical simulations, from which the cuspy behavior is 
deduced for the inner parts by means of extrapolation, 
are actually testing the halo shape down to 0(0. 1) kpc 

In our analysis we will take into account this large 
uncertainty in the inner core radius by varying r cut in 
the range [10~ 8 , 10" 1 ] kpc. 

C. Results for $ cosm ° 

The results of the calculations of the cosmological fac- 
tor $ cosmo for the Milky Way are shown in Fig. [3 for 
the four main profiles previously discussed and for a de- 
tector with angular resolution equal to 1° and 0.1°. A 
constant-density central region of radius r cut = 10~ 8 kpc 
has been used for the cuspy profiles. Since the value 
of r cu t used in Fig|2 somehow represents a lower bound 
on the acceptable values of this parameter, the values of 
•jcosmo shown in Fig. UJcan be taken as an upper bound 
on the cosmological factor, for any given halo profile and 
for the two representative acceptance angles. Clearly the 
non-cuspy profiles are not affected by the choice of r cu t- 

In the same figure, the values of $ cosmo for LMC and 
M31 are also shown. We see that these external galaxies 
can be resolved against the galactic signal in all cases, ex- 
cept for the case of LMC with an iso-core density profile. 
These two external galaxies can therefore be looked at as 
gamma-ray sources from DM annihilation (provided that 
the ensuing gamma-ray flux can be detected against the 
gamma-ray background). If a gamma-ray signal were de- 
tected, for instance from the galactic center, it should be 
correlated to a corresponding signal both from LMC and 
from M31. Since the "supcrsymmetric factor" is the same 
for all the sources, the relative strength of the gamma- 
ray fluxes from the galactic center, LMC and M31, could 
then be used to deduce information on the halo shape, 
since it depends only on the DM density profile. How- 
ever, this possibility is strongly limited by the fact that 
$ cosmo for LMC and M31 is much smaller than the one 
from the galactic center, as is clear from Fig. The en- 
suing fluxes from external galaxies will therefore be much 
smaller than the ones from the galactic center. 

The dependence of the cosmological factor on the cut- 
off radius of the inner core is shown in Fig. [3]for the Milky 
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FIG. 2: The lines denote the "cosmological factor" $ cosmo for 
the Milky Way, calculated for different dark matter profiles, 
for a solid angle AQ = 10 -3 sr (upper panel) and AQ — 10 _s 
sr (lower panel). The small boxes show a zoom at small angles 
toward the galactic center. A constant-density central region 
of radius r cu t = 10 -8 kpc has been used for the cuspy profiles. 
The points at ip ~ 81° and ip ~ 119° denote the values of 
$ cosmo for LMC and M31, respectively. From top to bottom 
the points refer to different halo profiles: Moore, NFW97, 
M04, iso-core in the upper panel; Moore, M04, NFW97, iso- 
core in the lower panel. 

Way and in Fig. 0]for LMC and M31. In these figures we 
plot the ratio $ cosmo (profile, r cut , i/')/$ cosmo (M99, r cut = 
lO- 8 kpc,-0 = 0) for the M99, NFW97, M04 and iso-core 
profiles and for r cut in the range discussed above . As 
expected, a cored distribution or a less cuspy profile than 
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FIG. 3: Relative strength of the line-of-sight integral with 
respect to different halo profiles and different inner core radii 
for the Milky Way. Numbers are normalized to the highest 
value of the $ cosmo given by a M99 profile with a physical 
cut-off radius of 10 _s kpc and at ip = 0. Left panel: solid 
angle AS7 = 10~ 3 sr. Right panel: solid angle AQ = 10 -5 sr. 

the M99 decreases the cosmological factor by a significant 
amount with respect to the most optimistic hypotheses of 
a M99 profile with an inner core radius r cu t = 1 CP 8 kpc. 
Figs. 03 and 0] quantify this effect. 

In the case of the Milky Way, the reduction factor at 
the galactic center can be sizeable: for instance, when a 
NFW97 profile with r cut = 0.1 kpc is used, the reduction 
is of the order of 4 • 10 -3 for a solid angle of observation 
AO = 10" 3 sr and 6 • 1CT 5 for AO = 10" 5 sr. In the 
case of the iso-cored distribution the reduction factor is 
as large as 1CT 4 for AO = lfT 3 and lfr 6 for Ail = lfT 5 . 

The same trend is observed for the external galaxies 
which we have considered, although the net effect is less 
prominent. In the case of M31, the reduction is at most 
a few xl0~ 2 for AO = 10 -3 and it can reach 10~ 3 for 
AO = 10~ 5 and the iso-core profile. For LMC, the re- 
duction is again of the order of 10 _2 -10 -3 , except for 
the iso-core profile and AO = 10~ 5 , for which it reaches 
values of the order of 10~ 5 . 

In the following, for definiteness we will refer to the 
most optimistic values of $ cosmo shown in Fig. [21 ob- 
tained for a M99 profile with a cut-off radius of 10 -8 kpc 
and to a NFW97 shape, with the same cut-off radius. Re- 
sults for different halo profiles or core parameters can be 
easily obtained by scaling the results according to Figs. 
and H 



FIG. 4: Relative strength of the line-of-sight integral with 
respect to different halo profiles and different inner core radii 
for M31 (upper panels) and the LMC (lower panels). Numbers 
are normalized to the highest value of the <j> cosmo given by a 
M99 profile with a physical cut-off radius of 10 -8 kpc and at 
tp = 0. Left panels: solid angle Af2 = 10~ 3 sr. Right panels: 
solid angle AO = 10~ 5 sr. 

D. Including Substructures 

In the CDM scenario, sub-halos that accrete into larger 
systems are tidally stripped of a fraction of their mass, 
originating debris streams Their dense central cores, 
however, survive the merging event and continue to orbit 
within the parent halo. High resolution N-body simula- 
tions [lj, |34| have indeed shown that DM halos host a 
population of sub-halos with a distribution function de- 
pending on the sub-halo mass and on the distance of the 
sub-halo from the halo center [3^. 

The effect of including sub-halos in the Milky Way and 
in the galaxies of the Local Group has been discussed in 
Refs. [llL I3H I37I . where different parameters for the 
sub-halo distribution, along with the existence of mass 
stripping and tidal heating, have been considered, and a 
minimum mass of 10 6 Mq was assumed for the sub-halos. 
The existence of such a sub-halo population leads to av- 
erage boost factors for expected rates which depend on 
the modeling of the sub-halo distribution and on the den- 
sity profile, and can range from few unities to more than 
10 4 . In no case, however, the field of view toward the 
Galactic Center is affected, since in that region the grav- 
itational strengthening reduces the probability of finding 
sub-halos. The total effect of the presence of sub-halos 
in external galaxies is limited to a factor 2-5. 
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A discussion on the minimum mass of sub-halos in 
our Galaxy can be found in Ref. [^3, where small scale 
clumps are considered, with masses down to 1O _8 M0 for 
a DM constituted by neutralinos. An average enhance- 
ment factor of 2-5 is found, depending on the profile, 
while the enhancement toward the Galactic Center is 
found to be of a factor - 0.3 (NFW97) to - 0.5 (M99). 

Hereafter we consider values for the "cosmological fac- 
tor" related to an unclumpy scenario. For a clumpy halo 
our results can be scaled according to the previous con- 
siderations. 



IV. THE "SUPERSYMMETRIC FACTOR" 

In our study we employ the MSSM supersymmetric 
extension of the Standard Model, which is defined as an 
effective theory at the electroweak scale. The scheme 
is defined in terms of a minimal number of parameters, 
only the ones which are necessary to shape the essen- 
tial properties of the theoretical structure of the MSSM 
and of its particle content. A number of assumptions are 
therefore imposed at the electroweak scale: a) all squark 
soft-mass parameters are degenerate: = m^; b) all 
slepton soft-mass parameters are degenerate: m ; ~. = mf 
c) all trilinear parameters vanish except those of the third 
family, which are defined in terms of a common dimcn- 
sionless parameter A: = = Arriq and Af = Ami. 
In addition, we employ also the standard relation at the 
electroweak scale between the U(l) and ST/ (2) gaugino 
mass parameters: Mi = (5/3) tan 2 Qw ~ 0.5 M2, which 
is the low energy consequence of an underlying unifica- 
tion condition for the gaugino masses at the GUT scale. 
In this class of gaugino-universal models, the neutralino 
mass has a lower bound of about 50 GeV. This limit is 
induced by the lower bound on the chargino mass deter- 



mined at LEP2 |39j: m x ± > 100 GeV. This is at variance 
with respect to effective MSSM schemes which do not 
posses gaugino-universality, where the neutralino mass 
can be as low as few GeV's (see for instance Refs. 
and references quoted therein). Gamma-ray detection 
from the annihilation of these light neutralinos has also 
been analized in Ref. jioj . in the context of SUGRA 
models where gaugino non universality is defined at the 
GUT scale. 

Due to the above mentioned assumptions, the super- 
symmetric parameter space of our scheme consists of the 
following independent parameters: M 2 , fi, tan/3, itia, 
rriq, mj and A. In the previous list of parameters fi de- 
notes the Higgs mixing mass parameter, itia is the mass 
of the CP-odd neutral Higgs boson and tan (3 = v t jvb is 



the ratio of the two Higgs v.e.v.'s that give mass to the 
top and bottom quarks. 

When we perform a numerical random scanning of 
the supersymmetric parameter space, we employ the fol- 
lowing ranges for the parameters: 1 < tan/3 < 50, 
100 GeV < M 2 < 6000 GeV, 100 GeV < m^mj < 
3000 GeV, sign(» = -1, 1, 90 GeV < m A < 1000 GeV, 
— 3 < A < 3. The range on both M 2 and [i extends up 
to 6 TeV in order to allow us to study also very heavy 
neutralinos, with a mass up to about 3 TeV. 

The parameters space of our effective MSSM is con- 
strained by many experimental bounds: accelerators 
data on supersymmetric and Higgs boson searches |4lj 
and on the invisible width of the Z boson, measure- 
ments of the branching ratio of the b — > s + 7 de- 
cay and of the upper bound on the branching ratio of 
B s — > /i + + fi~ , measurements of the muon anoma- 
lous magnetic moment a M = (g^ — 2)/2. The limits we 
use are: 2.18 • 10~ 4 < BR(b -> s + 7) < 4.28 • 10" 4 
0; BR(B S -> n+ + fj,~) < 7.5 • 10" 7 (95% C.L.) 0; 
-142 < Aa M • 10 11 < 474 (this 2a C.L. interval takes into 
account the recent evaluations of Refs. 0,0]). 

For the theoretical evaluation of BR(b — > s + 7) and 
BR(B S — > m + we have used the results of Ref. 0] 

and Ref. [47J, respectively, with inclusion of the QCD ra- 
diative corrections to the bottom-quark Yukawa coupling 
discussed in Ref. We notice that gluinos do not en- 
ter directly into our loop contributions to BR{b — > s + 7) 
and BR(B S — > /i + +/i~), since we assume flavor-diagonal 
sfermion mass matrices. Gluinos appear only in the 
QCD radiative corrections to the b Yukawa coupling: in 
this case M3 is taken at the standard unification value 
M 3 = M 2 a 3 (M z )/a 2 (M z ), where a 3 (M z ) and a 2 (M z ) 
are the SU(3) and SU(2) coupling constants evaluated at 
the scale M z . 

Another relevant observational constraint comes from 
Cosmology. The recent observations on the cosmic mi- 
crowave background from WMAP 0| , used in combina- 
tion with galaxy surveys, Lyman-a forest data and the 
Sloan Digital Sky Survey Collaboration results 0], are 
leading to a precise knowledge of the cosmological pa- 
rameters, and in particular of the amount of dark matter 
in the Universe. From the analysis of Ref. 0|, we can 
derive a restricted range for the relic density of a cold 
species like the neutralinos. The density parameter of 
cold dark matter is bounded at 2a level by the values: 
(Cl C DMh 2 ) min = 0.095 and (Sl CMf A 2 ) mH = 0.131. This 
is the range for CDM that we consider in the present 
paper. For supersymmetric models which provide val- 
ues of the neutralino relic abundance fl x h 2 smaller than 
the minimal value (£lcDMh 2 ) m i n , i.e. for models where 



the ncutralino represents a subdominant DM compo- 
nent, we accordingly rescale the value of the DM density: 

Px( r ) = &( r ) witn £ = ^x^ 2 /(^CI5M/i 2 )min- 

We recall that the relic abundance fl x h 2 is essentially 
given by fl x h 2 oc (a^v)^, where (<J ann v} in t is the 
thermal-average of the product of the neutralino anni- 
hilation cross section and velocity, integrated from the 
freeze-out temperature in the early Universe down to the 
present one. The analytical calculation of <7 a nn relies on 
the full set of available final states: fermion-antifcrmion 
pairs, gluon pairs, Higgs boson pairs, one Higgs boson 
and one gauge boson, pairs of gauge bosons |5l| . We 
have not included coannihilation |52( in our evaluation of 
the neutralino relic abundance, since in our effective su- 
persymmetric model a matching of the neutralino mass 
with other particle masses is usually accidental, and not 
induced by some intrinsic relationship among the differ- 
ent parameters of the supcrsymmetric model, like instead 
in a constrained SUGRA scheme, The inclusion of coan- 
nihilation would not change the main results of our anal- 
ysis, since it would only reflect in a limited reshuffle of a 
small fraction of the points of the scatter plots displayed 
in the next Sections. 

A. The annihilation cross section 

As already stated in Sec. ^ the gamma-ray flux pro- 
duced by ncutralino annihilation depends on the thermal 
average of the neutralino self-annihilation cross section 
(cann^) in the galactic halo at present time. The be- 
haviour of (cTunnv) /m 2 ,, which is a relevant quantity in 
the calculation of the gamma-ray flux, is shown in Fig. 
|S] as a function of the neutralino mass and for the ef- 
fective MSSM we are using. We remind that {cr ann v) 
in general is different from (crann^)int 

which is responsi- 
ble for the determination of the relic abundance. The 
two cross sections closely follow each other only for s- 
wave annihilation. An inverse proportionality between 
the gamma-ray signal and the relic abundance is there- 
fore usually a good approximation, although deviations 
are present. This effect is shown in the box-insert in Fig. 
□ 

Other key ingredients for the determination of the 
gamma-ray signal are the branching ratios of the anni- 
hilation cross section into the different final states. For 
neutralino lighter than 1 TeV the branching ratios were 
shown in Ref. pSflj . Fig. extends the behaviour of the 
branching ratios for neutralino masses higher that one 
TeV. We see that in this case the dominant channels are 
the two gauge bosons and the gauge boson+Higgs boson 
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FIG. 5: The thermally-averaged annihilation cross-section 
divided by the square of the neutralino mass m x as a func- 
tion of m x in the frame of the eMSSM. Crosses show the 
WMAP-preferred zone for a DM dominant neutralino. In the 
small box the annihilation cross section at the present epoch 
is shown as a function of the neutralino relic abundance. 

final states. 



B. The Photon Spectrum 

The diffuse photon spectrum from neutralino annihi- 
lation originates from the production of fcrmions, gauge 
bosons, Higgs bosons and gluons. Both gauge bosons 
and Higgs bosons eventually decay into fcrmions. The 
hadronization of quarks and gluons, in addition to ra- 
diative processes, can produce 7-rays. The main chan- 
nel of production of 7-rays goes through the production 
and subsequent decay of neutral pions. The contribu- 
tion to the 7-ray spectrum from production and decay 
of mesons other than pions (mostly 77, 7/ , charmed and 
bottom mesons) and of baryons is usually subdominant 
as compared to ir° decay and it has been neglected. Neu- 
tralino annihilation into lepton pairs can also produce 
7-rays from electromagnetic showering of the final state 
leptons. This process can be dominant for E 1 <, 100 
MeV, when the neutralino annihilation process has a siz- 
able branching ratio into lepton pairs. In the case of 
production of r leptons, their scmihadronic decays also 
produce neutral pions, which then further contribute to 
the gamma-ray flux. 

As discussed in Ref. Q, we have evaluated the 
gamma-ray fluxes originating from hadronization and ra- 
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FIG. 6: Branching ratios for high mass neutralino annihi- 
lation into fermions (upper left panel) and relative strength 
of annihilation into gauge bosons (upper right), Higgs bosons 
(lower left) and a gauge boson and a Higgs boson (lower right) 
with respect to the annihilation into fermions. 



FIG. 7: The photon spectrum from a m x — 1 TeV neutralino 
annihilation into: (a) leptons, (b) gauge bosons, (c) Higgs 
bosons decaying into r's and (d) Higgs bosons decaying into 
b's. For each curve a branching ratio of 100% in that channel 
has been considered. 



diativc processes by means of a Monte Carlo simulation 
with the PYTHIA package [§4|. In the present paper 
we extend that analysis by giving explicit fits to our nu- 
merical distributions which are valid for the energies of 
interest in the current analysis, i.e. for photon energies 
E > 10 GeV. When a flux is presented for energies below 
10 GeV, the numerical analysis has been used. 

The differential spectra of photons from DM annihila- 
tion have been parametrized as follows: 



dN. 
d.r 



(10) 



where x = E~ l /m x and i identify quarks, W, Z and glu- 
ons. The value of r\ is 2 for W, Z and top quark final 
states, and 1 otherwise. In the case of r leptons, the 
functional form for the differential number of photons is: 



dx 



x aT (b T x + c T x 2 + d T x 3 )e e 



(11) 



The values of the parameters of the fits are given in 
Tables [V] and IVII for the two representative values of 
m x = 500 GeV and m x = 1 TeV. 

In Fig. [3we show an example of photon spectra origi- 
nated by neutralino annihilation into different pure final 
states of a neutralino with m x = 1 TeV. We see that at 
lower energies the dominant contribution is given by the 
7-rays coming from the hadronization of quarks and glu- 



ons. The spectra coming from gauge bosons arc some- 
what harder, while the hardest ones are given by the 
r lepton. In the case of Higgs bosons, the spectra are 
mainly driven by the type of particle in which the Higgs 
bosons decays, and are somewhat softer. 



C. Resuls on <E> SUSY 

The composition of the information on the neutralino 
annihilation cross-section and branching ratios with the 
informations coming from the differential spectra of pho- 
tons from the annihilation of neutralinos in pure final 
states provides the prediction of what we have called the 
"supersymmetric factor" $ SUSY i n the 7-ray flux compu- 
tation. Fig. shows <I> SUSY defined as the integral of the 
gamma-ray flux of Eq. @ above a set of sample thresh- 
old energies: 1, 10, 50 and 100 GeV. When the threshold 
energy is small, $ SUSY is rougly inversely proportional 
to the neutralino mass. Since the neutralinos annihi- 
late almost at rest, when the threshold energy increases 
(jjSUSY s ig n i ncan tly drops because the highest available 
photon energy is E ~ m x for any given neutralino mass. 
The highest value for $ SUSY is of the order of 10~ 8 cm 4 
kpc -1 s _1 GeV" 2 sr" 1 when the threshold energy is 1 
GeV. At masses larger than about 500 GeV for any given 
threshold energy the values of $ SUSY all lie inside a band 
which span no more than one order of magnitude. This 
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TABLE V: Fitted parameters of Eq. 1101 for the annihilation 
of neutralinos into quarks and gauge bosons, calculated for 
m x — 500 GeV and m x = 1 TeV. Fits obtained with these 
parameters are valid down to E = 10 GeV. 





m x = 500 GeV 


m x = 1 TeV 


a T 


-1.34 


-1.31 


br 


6.27 


6.94 


Ct 


0.89 


-4.93 


d T 


-4.90 


-0.51 


e T 


-5.10 


-4.53 



TABLE VI: Parameters of Eq. (I I I ti for the annihilation of 
neutralinos into r leptons, calculated for m x = 500 GeV and 
m x = 1 TeV. Fits obtained with these parameters are valid 
down to E = 10 GeV. 



makes the predictions on the gamma-ray fluxes for large 
neutralino masses quite predictive: the possible variation 
due to the different supersymmctric models is confined to 
a relatively restricted range, much smaller than for the 
case of lighter neutralinos. 

The information on the factor <J> SUSY is detailed in Ta- 
bles I VIII IVIIII and IIXI where we give the number of pho- 
tons produced in each pure final state for different thrcsh- 
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FIG. 8: The "supersymmetric factor" $ as a function of 
the neutralino mass. Different colors show different threshold 
energies above which the energy spectra have been integrated. 

old energies. This information may be used to make pre- 
dictions for the gamma-ray fluxes also for DM candidates 
other than the neutralino. 

The results of this Section and of Sec. IIIII will be used 
in the next Sections to predict the photon fluxes from 
the galactic center and from our representative external 
galaxies. 

V. PREDICTION AND DETECTABILITY OF 
PHOTON FLUXES 

In this Section we will show the results on the pre- 
diction of photon fluxes from neutralino annihilation in 
our Galaxy and in some selected external galaxies. We 
will therefore study the detectability of such signals with 
ground-based Cerenkov telescopes and next generation 
satellite-borne experiments. 

A. Predicted Photon Fluxes from Neutralino 
Annihilation 

In the previous Sections we have computed the "cos- 
mological factor" <j> cosmo (see Fig. EJ and the "supersym- 
metric factor" $ SUSY (see Fig. |HJ. We are now ready 
to predict the gamma-ray fluxes fron neutralino anni- 
hilation in the effective MSSM. Results are reported in 
Figs. El and ^3 where we show the expected fluxes of 
7-rays with energies above 50 GeV and 100 GeV from 
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TABLE VII: Integrated number of photons above a given en- 
ergy E from the annihilation of neutralinos with masses 500 
GeV, 800 GeV and 1 TeV, for different channels of annihila- 
tion. 



the galactic center and M31 and for a detector aperture 
of Af2 = 10~ 5 sr. Fig. [5| refers to the galactic center 
for a Milky Way with a NFW97 density profile, while 
Fig. is calculated for M31 with a M99 density pro- 
file. The spread of points is given by the different SUSY 
parameters corresponding to each point. 

In the case of the flux from the galactic center with 
a NFW97 profile and a typical threshold energy of 50 
GeV, we predict a maximal gamma-ray flux of the or- 
der of 10~ 12 cm" 2 s" 1 for neutralinos lighter than 200 
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TABLE VIII: Integrated number of photons above a given 
energy E from the annihilation of neutralinos with masses 500 
GeV, 800 GeV and 1 TeV, for different channels of annihila- 
tion. 



GeV, while heavier neutralinos can provide a maximal 
flux of the order of a few 10" 13 cm" 2 s" 1 . In the case 
of a M99 density profile toward the galactic center, the 
fluxes are increased by a factor of about 160, as can be 
deduced from Fig. [31 In this case the maximal fluxes 
can reach the level of 10~ 10 cm" 2 s" 1 . If the detector 
threshold energy is increased to 100 GeV the gamma- 
ray fluxes are one order of magnitude smaller. Finally, 
as a consequence of the previously discussed property of 
$ su , we see that for neutralino masses heavier than 
about 500 GeV the supersymmetric models we are con- 
sidering provide gamma-ray fluxes inside a band with a 
lower limit of a few 10" 14 cm" 2 s" 1 , for a NFW97 pro- 
file. Obviously, if we enlarge the allowed intervals for 
the MSSM parameters (our definitions are given in Sec. 
II V|) . lower gamma-ray fluxes can be obtained also for 
heavy neutralinos. However, if we consider natural mass 
scales for the supersymmetric model, which means that 
we should not increase the scale of the mass parameters 
of the model much over the TeV scale, Fig. [21 shows the 
level of the lower limit on the gamma-ray flux for heavy 
neutralinos. 

Also the Andromeda Galaxy can provide gamma-ray 
fluxes of the order of 10 _12 -10~ 13 cm -2 s" 1 inside a 
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TABLE IX: Integrated number of photons above a given 
energy E from the annihilation of 1 TeV neutralino into a 
sample state of Higgs bosons, with subsequent decay into 6 
quarks or tau leptons. A mass of 120 GeV has been assumed 
for the light Higgs, while a mass of 500 GeV has been taken 
for the charged, heavy and pseudoscalar Higgses. 



solid angle of Af2 = 10~ 5 sr, but only for a M99 den- 
sity profile. These values therefore represent the maximal 
fluxes which can be produced by neutralino annihilation 
in M31. We remind that although the galactic center is 
much brighter for the same density profile, M31 can be 
resolved over the galactic gamma-ray signal due to its 
location at ip = 119°, as is shown in Fig. [3 

In the following we will compare our expected fluxes 
with the sensitivity curves of foreseeable experiments. 
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FIG. 9: Integrated gamma-ray fluxes from neutralino anni- 
hilation at the galactic center, for a NFW97 density profile 
and inside a solid angle AQ = 10~ J sr. Two representative 
threshold energies have been assumed: 50 GeV (upper panel) 
and 100 GeV (lower panel). 



B. Detectability of Photon Fluxes from Neutralino 
Annihilation 



We have considered two platforms of observations of 
7-rays from neutralino annihilation, corresponding to a 
Cerenkov apparatus with the characteristics of VERITAS 
1 1 and to a satellite-borne experiment similar to GLAST 
4J. The detectability of the diffuse flux from DM an- 
nihilation is computed by comparing the number n 7 of 
expected 7 events with the fluctuations of background 
events ribkg- To this purpose we define the following ra- 



tio a given by: 



(12) 



/ Af(E, e)[d(j)^ M /dEdVL]dEdVL 



All s (E, 6)[d<t> bkg /dEdn]dEdn 

where T$ defines the effective observation time and <f>bkg 
is the background flux. For a Cerenkov apparatus, for in- 
stance, it is defined as the time during which the source is 
seen with zenith angle 9 < 60°. The quantity £aq = 0.7 
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together more Cerenkov telescopes, the latter is intrinsi- 
cally limited by the size of the satellite and cannot be 
much greater than the fiducial value quoted here. Fi- 
nally we have assumed an angular resolution of 0.1° for 
both instruments, and a total effective pointing time of 
20 days for the Cerenkov telescope and 30 days for the ex- 
periment on satellite. An identification efficiency e must 
be taken into account, which is one of the most important 
factors which have to be studied in order to reduce the 
physical background level. A Cerenkov apparatus has 
a typical identification efficiency for electromagnetic in- 
duced (primary 7 or electrons) showers e c . m . ~ 99% and 
for hadronic showers ehad ~ 99%. This means that only 1 
hadronic shower out of 100 is misidentified as an electro- 
magnetic shower. Unfortunately, this method cannot dis- 
tinguish between primary photons and electrons, which 
therefore represent an irreducible background for ground- 
based detectors. As far as a satellite-borne experiment 
is concerned, an identification efficiency for charged par- 
ticles of e c har g ed ~ 99.997% can be assumed, while for 
photons it lowers to e ne ujgal ~ 90% due to the backsplash 
of high energy photons [55| ■ 

We have considered the following values for the back- 
ground levels. For the proton background we use |56j : 



1A9E 



-2.74 



P 



dQdE cm 2 s sr GeV ' 

while for the electron background j^: 

-2 ci-3.3 e 



= 6.9 x lQ'^E 



(13) 



(14) 



dQ.dE cm 2 s sr GeV 

and finally for the Galactic photon emission, as extrapo 
lated by EGRET data at lower energies, we employ |58j 



FIG. 10: Integrated gamma-ray fluxes from neutralino anni- 
hilation in M31, for a M99 density profile and inside a solid 
angle Afi = 10~ J sr. Two representative threshold ener- 
gies have been assumed: 50 GeV (upper panel) and 100 GeV 
(lower panel). 



is the fraction of signal events within the optimal solid 
angle Afl corresponding to the angular resolution of the 
instrument. The effective detection areas A for electro- 
magnetic and hadronic induced showers are defined as the 
detection efficiency times the geometrical detection area. 
For the case of a Cerenkov apparatus we have assumed a 
conservative effective area A cS = 4 x 10 8 cm 2 , while for a 
satellite experiment we have considered A = 10 4 cm 2 . 
Both values have been assumed independent from E and 
9. Note that while the former can be increased by adding 



i.gal-7 
^diffuse 



N (l,b) 10" 6 £ Q 



7 



(15) 



dfldE --"v."/" 7 cm 2 s sr GeV ' 

with a set to —2.7 in all the considered energy range, in 
lack of data for energies higher than tens of GeV. The 
normalization factor Nq depends only on the interstellar 
matter distribution, and is modeled as |5£ 



N (l,b) = 



85.5 



1 + (Z/35) 2 J I + (6/(1.1 + |Z| 0.022)) 2 



+ 0.5 
(16) 



for |/| > 30° and 
N (l,b) 



85.5 



^l + (Z/35) 2 ^/l + (6/1.8) a 



0.5 (17) 



for | Z | < 30°, where the longitude I and the latitude b 
are assumed to vary in the intervals —180° < I < 180° 
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and —90° < b < 90°, respectively. Finally, for the diffuse 
extragalactic 7 emission, as extrapolated from EGRET 
data at lower energies [59j, we use: 



N. Forrengo. L. Fieri. S. Scepe (2CC4I 

Galactic Centre 



m extra— 7 
'^diffuse 

dndE 



1.38 x 1(T D £; 



-6j?-2.1 



7 



cm 2 s sr GeV 



(18) 



If a galactic origin of high galactic latitude 7 emission is 
considered, then this last estimate should be increased 
by about 60% 0. 

Fig. ^Jshows the 5 a sensitivity curves for the experi- 
mental apparata discussed above. Due to the different 7 
backgrounds, the curves are slightly different in the direc- 
tion of the galactic center or toward the M31 galaxy. Also 
plotted for reference is the expected integrated 7-ray flux 
for a SUSY model with m x = 1 TeV, 50% branching ratio 
of annihilation into W bosons and 50% into Higgs bosons 
(following the results of Fig. Elfor the branching ratios of 
high mass ncutralinos), and an annihilation cross-section 
of 2 x 10~ 26 cm 3 s~ : which refers to the most optimistic 
values of Figs. l9landHHI Due to our discussion in the pre- 
vious Section on the properties of $ SUSY ; one could then 
consider the curve of 7-ray flux from neutralino annihi- 
lation which we show in Fig. II II as the highest spectrum 
of a range of curves given by the spread of points in Figs. 
EDandUni 

From Fig. 1111 and our previous discussion on the cos- 
mological and supersymmetric factor it therefore arises 
that signals from extragalactic objects could hardly be 
detected. The gamma-ray spectrum calculated for a M99 
profile is two orders of magnitude smaller than the ex- 
pected sensitivities we estimate for detectors like GLAST 
and about one order of magnitude smaller than the esti- 
mated sensitivity of VERITAS. We also notice that the 
most optimistic prediction for the flux we are showing in 
Fig. ^]is at the level of the extrapolated background, a 
fact which by itself would make problematic the obser- 
vation of a signal from M31. Only in the very optimistic 
case of a clumpy M99 matter density, the expected signal 
would exceed the extrapolated background, but it would 
nevertheless remain unaccessible. 

In the case of a signal from the galactic center, a den- 
sity profile as cuspy as M99 (or the adiab-NFW) could 
be resolved by both a satellite detector like GLAST and 
a Cerenkov telescopes with the characteristics of VERI- 
TAS. In the case of a NFW97 profile, a potential signal 
would not be accessible. Therefore, in the case of the sig- 
nal from the galactic center a density profile harder than 
NFW97 is required in order to have a signal accessible to 
GLAST-like and VERITAS-like detectors. 
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FIG. 11: Study of the sensitivity of an ACT detector and 
a satellite borne experiment to photon fluxes from a TeV 
neutralino annihilation. Solid lines denote the 5a sensitiv- 
ity curves for satellite and Cerenkov detectors. These curves 
have been calculated according to the prescriptions given in 
the text. The flux expected from the GC with a NFW97 and 
a M99 profile are shown in the upper panel. The flux from 
M31 with a M99 profile is shown in the lower panel. Pho- 
ton fluxes are given for Afi = 1CF 5 sr, which is the typical 
detector acceptance. 



C. Comparison with Recent Data 

Recent experimental data taken from CANGAROO-II 
in the direction of the galactic center, show that the 
spectral shape of photons from the GC is in excess of the 
extrapolated background from standard processes. Fig. 
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IHlshows the CANGAROO-II data in the right panel, and 
the EGRET data [jj at lower energies in the left panels. 
We have superimposed to the data the 7-ray background 
used in our previous analysis, as well as the predicted 
7-ray spectra from high mass neutralino annihilation, for 
the NFW97 and the M99 profiles. These spectra have 
been normalized within a solid angle coherent with the 
observations. We can see that not even a M99 profile can 
reproduce the observed data, as already observed in Ref. 
[61|. Fig. EJrcproduces the same information of Fig. ^1 
but the "cosmological factor" has been enhanced by a fac- 
tor 2.5 (equivalently, one could think to an enhancement 
in the "supersymmetric factor" , but this is not possible in 
the effective MSSM, neither in more constrained minimal 
SUGRA models which usually provide annihilation cross 
sections smaller than the effective MSSM). We can see 
that, when appropriately boosted, the signals from anni- 
hilation of ncutralinos with mass higher than 1 TcV have 
the property of matching the observed CANGAROO-II 
data and not being in conflict with the EGRET data. 

On the other hand, Fig. H2lshows that it is not possible 
to explain at the same time both the EGRET excess in 
the 1-20 GeV energy range and the CANGAROO-II flux 
at energies above 250 GeV with the spectral shape of a 
gamma-ray flux from neutralino annihilation. While the 
EGRET spectrum can be well explained by a light neu- 
tralino in a non-universal gaugino model pj , with niy ~ 
30-40 GeV, or by a neutralino of about 50-60 GeV 0| 
in the effective MSSM, the CANGAROO-II data require 
much heavier neutralinos in order to produce photons in 
the hundred of GeV range: in this case, however, the en- 
suing gamma-ray spectra are too low in the 1-10 GeV 
range and cannot reproduce the EGRET data together 
with the CANGAROO-II ones. 

We complete this Section by applying our method to 
M87 and comparing our results with the measurements 
available for that galaxy, which show a possible indication 
of a 7-ray excess. This is shown in Fig. 1141 where one 
can see that our predictions are well below the flux mea- 
sured by HEGRA jf^] , even if a M99 profile is assumed. 
Not even a clumpy distribution, which could enhance the 
predicted fluxes by at most a factor of 5, would allow us 
to explain the HEGRA excess by means of neutralino 
annihilations in the effective MSSM. 

VI. CONCLUSIONS 

We have discussed the gamma-ray signal from dark 
matter annihilation in our Galaxy and in external ob- 
jects, namely the Large Magellanic Cloud, the An- 



N. Fernerigo, L. Fieri, 5. Scope (2004) 




10 3 



Energy (GeV) 

FIG. 12: Differential spectrum of the photon flux expected 
from neutralino annihilation in the galactic center. A 50% 
branching ratio into W pairs and 50% into b quarks has been 
assumed. Solid lines represent the calculation for a M99 
profile for different neutralino masses, while dashed-dotted 
lines show the same spectra assuming a NFW97 profile. Dot- 
ted lines shows the extrapolated 7-ray "conventional" back- 
ground. Open circles (left panel) show the EGRET results 
on photon flux from the galactic center, while filled circles 
(right panel) show the recent data at higher energies from 
CANGAROO-II. Photon fluxes are given for the correspond- 
ing typical detector acceptance, that is for Afi = 10 -3 sr in 
the left panel and for AST2 = 5 ■ 10 -5 sr in the right panel. 

dromeda Galaxy (M31) and M87. The aim of our pa- 
per was to derive consistent predictions for the fluxes 
in a specific realization of supersymmetry, the effective 
MSSM, and to compare the predictions with the capa- 
bilities of new-generation satellite-borne experiments, 
like GLAST, and ground-based Cerenkov telescopes, for 
which we have used, for definiteness, the characteristics 
of the VERITAS telescope. 

Our results show that only the signal from neutralino 
annihilation at the galactic center could be accessible 
to both satellite-borne experiments and to ACTs, even 
though this requires very steep dark matter density pro- 
files toward the galactic center. A profile steeper than 
NFW97 is required in order to provide signals which can 
reach detectable levels. In the case of signals coming from 
external galaxies, even though the extragalactic signal is 
larger than the galactic contribution from neutralino an- 
nihilation, nevertheless the absolute level of the flux is 
too low to allow detection with the experimental tech- 
niques currently under development. 
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FIG. 13: The same as in Fig. [TJ]for a M99 profile multiplied 
by a factor 2.5 (dashed lines). 
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FIG. 15: The same as in Fig |12l including the data from 
HESS |i3 (see Note Added at the end of the paper). Photon 
fluxes are shown for neutralino masses up to 20 TeV and for 
an M99 density profile. 

malization of the flux requires a boost factor of about 2.5 
over the flux obtained with a Moore et al. profile: this 
seems hard to obtain even in the presence of clumps. We 
also showed that the agreement with the CANGAROO- 
II data which is obtained with these boosted fluxes is not 
in contrast with the lower-energy EGRET data from the 
galactic center. In addition we showed that the spectral 
features of such fluxes cannot explain at the same time 
both the CANGAROO-II and EGRET excess by invok- 
ing a very heavy neutralino. Finally, we compared our 
predictions for the signal from M87 with the HEGRA 
data and found that the predicted fluxes from neutralino 
annihilation are too low to explain the HEGRA result. 



FIG. 14: Integrated photon flux as expected from a TeV 
neutralino annihilation in the M87 galaxy. Photon fluxes are 
given for Af2 = 10 -5 sr, which is the typical detector accep- 
tance. Also shown on the figure the upper limit determined 
by WHIPPLE 1(3 and the measurement from HEGRA 

We have also compared our theoretical predictions 
with the recent CANGAROO-II data from the galactic 
center and with the HEGRA data from M87. In both 
cases an indication of a gamma-ray excess is present. In 
the case of the CANGAROO-II data, the spectral shape 
is well reproduced by a gamma-ray flux from annihila- 
tion of neutralinos somewhat heavier than about 1 TeV, 
in agreement with Ref. [6l| . However the overall nor- 



VII. NOTE ADDED 

The HESS Cerenkov telescope j6^| has recently pub- 
lished new data on gamma rays from the galactic cen- 
ter. The measured flux and spectrum differ substan- 
tially from previous results, in particular those reported 
by the CANGAROO collaboration, exhibiting a much 
harder power-law energy spectrum, with spectral index 
of about —2.2. According to our analysis, these data, 
if interpreted in terms of neutralino annihilation, would 
require a neutralino mass in the range 10 TeV < m x < 
20 TeV and an M99 profile for the DM distribution, as 
shown in Fig |15l 
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